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Abstract. A computer-algebra aided method is carried out, for determining geometric 
objects associated to differential operators that satisfy the elliptic ansatz. This results in 
^ ' examples of Lame curves with double reduction and in the explicit reduction of the theta 

function of a Halphen curve. 
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1 Introduction 

The classical theory of reduction, initiated by Weierstrass, has found modern applications to 
^ ■ the fields of Integrable Systems and Number Theory, to name but two. In this short note 

we only address specific cases, providing minimal historical references, listing the steps that 
we devised, and exhibiting some new explicit formulas. A full-length discussion of original 
motivation, theoretical advancements and modern applications would take more than one book 
to present fairly, and again, we choose to provide one (two-part) reference only, which is recent 
and captures our point of view [31 0]. 

Our point of departure is rooted in the classical theory of Ordinary Differential Equations 
(ODEs). At a time when activity in the study of elliptic functions was most intense, Halphen, 
Hermite and Lame (among many others) obtained deep results in the spectral theory of lin- 
ear differential operators with elliptic coefficients. Using differential algebra, Burchnall and 
Chaundy described the spectrum (by which we mean the joint spectrum of the commuting 
operators) of those operators that are now called algebro-geometric, and some non-linear Par- 
tial Differential Equations (PDEs) satisfied by their coefficients under isospectral deformations 
along 'time' flows, t\, . . . , t g , where g is the genus of the spectral ring. We seek algorithms that, 
starting with an Ordinary Differential Operator (ODO) with elliptic coefficients, produce an 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at http://www.emis.de/journals/SIGMA/kuznctsov.html 

1 Wishes to express thanks to the MISGAM project for support of his research visit to SISSA (Trieste) in 2006. 
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algebraic equation for an affine plane model of the spectral curve (Section [2|). This problem 
turns out to merge with Weierstrass' question [IB], Appendix C: 

Give a condition that an algebraic function f(x) must satisfy if among the integrals 



where F is a rational function of x and f(x), there exist some that can be transformed into 
elliptic integrals. 

This means that the curve X whose function field is C(x, f{x)) has a holomorphic differential 
which is elliptic, or equivalently, Jac X contains an elliptic curve. Weierstrass gave a condition in 
terms of the period matrix of X. Burchnall-Chaundy spectral curves of operators with elliptic 
coefficients automatically come with this reduction property, because the differential-algebra 
theory guarantees that all the coefficients of the spectral ring are doubly-periodic with respect 
to the same lattice; the PDEs mentioned above have solutions which are called 'elliptic solitons'. 
A more special condition is that there be a 'multiple' reduction, namely another periodicity 
with respect to a (possibly non-isomorphic) lattice, in which case a curve of genus three has 
Jacobian isogenous to the product of three elliptic curves, and this is what we found in [9]. 
Another difficult question is whether the Jacobian is actually complex-isomorphic (without 
the principal polarization, since Jacobians are indecomposable) to the product of three ellip- 
tic curves. This is the case for the Klein curve [2]; we do not know if it is for the Halphen 
curve in [9j. In genus two, the situation is only apparently simpler: the question of complex- 
isomorphism to a product has only recently been settled [6J, and that, in an analytic but not 
in an algebraic language; the question of loci of genus-2 curves which cover (n : 1) an ellip- 
tic curve remains largely open (despite recent progress, much of it based on computer algeb- 



In Section [3] of this paper, we provide progress towards such questions for the case of Lame 
curves, which are hyper elliptic: the curves are known, at least for small genus, but we re- 
port a double reduction (equianharmonic case) that we have not seen elsewhere; in Section [H 
we treat the genus 3 non-hyperelliptic case. Our contribution to a method consists of two 
parts. First, we provide algebraic maps between the higher-genus and the genus-one curves, by 
giving a search algorithm that identifies genus- 1 integrals, the very question posed by Weier- 
strass. Second, using classical reduction theory as revisited by Martens [E], we write the theta 
function for a case of reduction in terms of genus- 1 theta functions (or the associate sigma 
functions). We connect the two issues, by calculating explicitly the period matrix for some 
curve that admits reduction, given algebraically. It is a rare event, that a period matrix can 
be calculated from an algebraic equation, since the two in general are transcendental functions 
of each other. The method consists in the traditional one [2] [26"1 [27] of choosing a suitable 
homology basis and calculating the covering action on it. We believe our result on the period 
matrix of the Halphen curve and its reduction to be new. We also address the question of 
effectivization of the solution of the KP equation, or hierarchy; for this to be complete, how- 
ever, we need not only to reduce the theta function, but also the integrals of differentials of the 
second kind. 

We dedicate this small note to the memory of Vadim Kuznetsov. It was one of Vadim's many 
projects to find a systematic answer to the problem of giving a canonical form for differential 
equations that have a basis of solutions polynomial in elliptic functions, generalizing the theory 
of elliptical harmonics given in the last chapter of [36], and which we have implemented in 
Section [3] to study the Lame curves. In fact, one among us (VZE) co-authored one of the last 
papers of Vadim [37J, devoted to a problem that arises in physics and was first identified by 
him; some of the ideas that we describe below were to be used in the further development of 
that subject. The loss of Vadim as friend and mathematician has affected us deeply. Personal 




ra [29]). 
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memories and testimonials are collected in other parts of this issue; here we only say that it was 
our great privilege to know him. 

2 Elliptic covers 

Curves that admit a finite-to-one map to an elliptic curve ('Elliptic Covers') are special, and 
are a topic of extensive study in algebraic geometry and number theory. They have come to the 
fore in the theory of integrable equations, in particular their algebro- geometric solutions, for the 
reason that a yet more special class of such curves provides solutions that are doubly-periodic 
in one variable ('Elliptic Solitons'). 

Our point of view on the latter class of curves is based on the Burchnall-Chaundy theory, via 
which, we connect the problem to the theory of ODOs with elliptic coefficients, much studied 
in the nineteenth century (Lame, Hermite, Halphen, e.g.). It is only when such an operator 
has a large commutator that this theory connects with that of integrable equations. While the 
nineteenth-century point of view could not have phrased the special property in this way, the 
related question that was addressed was, when does the series expansion of solutions, in some 
suitable local coordinate on the curve (cf. Remark 1 in Subsection 3.1) terminate? The early 
work turned up important, but somehow ad hoc, properties that characterize algebro-geometric 
operators. 

2.1 Challenge I 

Detect (larger) classes of algebro-geometric operators. 

We do not address this issue here. A CAS-based project would be to test whether the operator 

r ^3 4 2 Q 2n(n-3)(4n + 3) , 
L = d — -n 2 pd + — i ^ V 

has centralizer larger than C[L]: we ask the question since Halphen |14j noticed that the equation 
L = has solutions that can be expressed in terms of elliptic functions. A more far-reaching 
project would construct elliptic-coefficient, finite-gap operators starting with curves that were 
shown to be elliptic covers, choosing a local coordinate whose suitable power is a function pulled 
back from the elliptic curve, to give rise to the commutative ring of ODOs, according to the 
Burchnall-Chaundy-Krichever inverse spectral theory. Curves of genus 3 that are known to be 
elliptic 3 : 1 and 2 : 1 covers, respectively, are Klein's and Kowalevski's quartic curves: 

zx 3 + xy 3 + yz 3 = 0, 

(z 2 — f2) 2 = 4xy(ax + by)(cx + aly) 

in projective coordinates, with fi a form of degree 2 in cc and y (such that the quartic is non- 
singular). Notably, Klein's curve (the only curve of genus 3 that has maximum number of 
automorphisms, 168) has Jacobian which is isomorphic to the product of three (isomorphic) 
elliptic curves [21 [231 125] ; Kowalevski, as part of her thesis, classified the (non-hyperelliptic) 
curves of genus 3 that cover 2 : 1 an elliptic curve (Klein's curve being a special case). 

2.2 Challenge II 

Find an equation for the spectral curve. 

This can be found once an algebro-geometric ODO is given. While the Burchnall-Chaundy 
theory allows one in principle to write the equation of the curve as a differential resultant |24j . 
this would be unwieldy for all but the simplest cases: first, one would have to solve the ODEs 
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of commutation [L, B] = for an unknown operator B of order co-prime with L (the simplest 
case of a two-generator centralizer) ; then, once B is known, calculate the determinant of the 
(m + n) x (m + n) matrix that detects the common eigenfunctions of L — A and B — fx [5j. 

Instead, it is much more efficient to adapt Hermite's (or Halphen's) ansatz; in the case of 
Lame's equation, the case of solutions 'written in finite form' is treated in [36]; more generally, 
assume that an eigenfunction can be written in terms of: 

¥(x;or) = e kx ^2aj(x,iJ,,k)--—<f>(x;a) 

3=0 X 

solving a given spectral problem: L^(x) = z^f(x), where: 
a{a)a{x) 

and a is a complex number viewed as a point of the elliptic curve v 2 = 4/i 3 — Q2[i — 53. 

By expanding p(x) and 4>(x, a) near x = and comparing coefficients, the aj can be written 
in terms of (k, z) and taking the resultant of the compatibility conditions gives an equation of 
the spectral curve in the (k, z)-plane [7]. 

2.3 Challenge III 

Express the eigenfunctions in terms of the theta function of the curve. 

2.4 Challenge IV 

Turn on the KP (isospectral-)time deformations ti and express the time-dependence of the 
coefficients of the operators (in particular, one such coefficient is an exact solution of the KP 
equation) . The expression of a KP solution in terms of elliptic functions is part of the programme 
sometimes referred to as 'effectivization'. 

Our approach to such 'challenges' is two-fold: on one hand we use the explicit expression 
of the higher-genus sigma functions that were classically proposed by Klein and Baker, but 
revived, substantively generalized, and brought into usable form at first by one of the authors 
with collaborators (cf. e.g. [5j). On the other hand, we design computer algebra routines which 
allow us to read, for example, coefficients of multi-variable Taylor expansions far enough to 
obtain the equations for the curves and the KP solutions for small genus. 

In this paper we present classes of examples, originating with classical ansatzes of Lame, 
Hermite and Halphen (specifically applied to elliptic solitons by Krichever [17]) as well as some 
detail of our general strategy. 



3 Hyperelliptic case 

The Lame equation has been vastly applied and vastly studied: we refer to [18] for general 
information, and to the classical treatise [36] for calculations that we shall need (for richer 
sources please consult pj~], and p2] for different perspectives). Our point of view is that the 
potentials n{n + l)p(£) are finite-gap when n is an integer (an adaptation of Ince's theorem to 
the complex-valued case) where the equation is written in Weierstrass form: 

f01-n(n+l)p(ZMO=MO (1) 
(for the rational form see [18] or [36J). 
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3.1 The spectral curve 

The Lame curves are hyperelliptic since one of the commuting differential operators has order 2, 
and will have genus n for the operator in ([1]). There is no 'closed form' for general n, but the 
equations have been found by several authors and methods, cf. specifically [31 d] and references 
therein (work by Eilbeck et al, by Enol'skii and N.A. Kostov, by A. Treibich and by J.-L. Verdier 
is there referenced), [13 |, 1181 l2"Tj 134] . In particular, [18] adopts a method of inserting ansatzes in 
the equation and comparing expansions similar to ours, and tabulates the equations for 1 < g < 8 
(in principle, a recursive calculation will produce them for any genus), so we do not present our 
table (1 < g < 10), but use g = 3 to exemplify our treatment of reduction in Subsection 3.2 
below. 

Remark 1. To motivate the project that Vadim Kuznetsov had wanted to design, as recalled in 
the introduction, we mention in sketch the steps of our strategy to find the spectral curves. The 
key ansatz is the finite expansion of the eigenfunctions; we used a formula from [36] in the case 
of odd n and the (equivalent) table in [1] for even n just for bookkeeping purposes; respectively, 

(n-3)/2 n/2 

w(0 = Vp(0-eWp(0-e2 £ Mp(£)-e 3 r /2 - r -\ w(0 = £ A r (-p(0) r . 

r=0 r=0 

It is then a matter of finding the coefficients b r , A r . For this, we insert w in (1) and expand in 
powers of p(£) or (p(£) — e^) 1 / 2 , depending on the various cases; we also substitute for the in 
terms of the modular functions g2 and g^. By comparing coefficients in the expansions we obtain 

a set of equations for the b r or A r , and we solve the compatibility condition of these equations. 

3 

The Lame curve f s (z) Y\ fi( z ) is the product of 4 factors, one f s (z) symmetric in the e^, i.e. just 

i=l 

depending on gi and 53, and three fi(z) which depend on one of the e, and $2 an d 93, which 
we found in turn. At the end, to be sure, for a given root z of one of these factors, we can 
determine the b r , resp. A r up to a common factor, and write the eigenfunction as a polynomial 
in p(£). The finiteness condition we used in expanding the eigenfunctions (a sort of Halphen's 
ansatz, cf. Section 1) characterizes finite-gap operators, in Burchnall-Chaundy terms, the ones 
that have large centralizers. It is worth quoting a speculation [28] p. 34 which regards finiteness 
in the local parameter z . 'We do not know an altogether satisfying description of the desired 
class C( n ); roughly speaking, it consists of the operators whose formal Baker functions converge 
for large z.' 

3.2 Cover equations and reduction 

Even though the Lame curves are elliptic covers a priori, it is not easy to find the elliptic 
reduction, which will correspond to a holomorphic differential of smallest order at the point 00 
in (w, z)-coordinates^]; by running through a linear-algebra elimination for a basis of holomorphic 
differentials, expressed in terms of rational functions of (w, z) (here again expansions provide 
the initial guess, lest the programs exhaust computational power), we exhibited the elliptic 
differential for all our curves. Note that this gives the degree of the cover, not known a priori. 
It is worth pointing out that more than n covers might be found (for genus n): indeed, according 
to a classical result (Bolza, Poincare, cf. [15] where a modern proof is given), if an abelian surface 
contains more than one elliptic subgroup, then it contains infinitely many. Note, however, that 
it can only contain a finite number for any given degree (the intersection number with a principal 
polarization) because the Neron-Severi group is finitely generated [15]. This degree translates 

2 A different and very efficient method was developed by Takemura in the beautiful series [311 32, 33, 34] where 
elliptic-hyperelliptic formulas are found by comparing the monodromies on the two curves. 
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into the degree of the elliptic cover. We are not claiming to have found all of them, or that those 
which we found be minimal. For example, for genus n = 2 (where the reduction is known since 
the nineteenth century), finding f s = z 2 — 3g 2 , fi = z + 3ej, in expanded form the curve is: 

w 2 = -(3(72 - z 2 ) (27(73 - 4z 3 + 9zg 2 ). 
The holomorphic differential is 

dp dz 

— = 3^— , 

p w 

and a second reduction (predicted by a classical theorem of Picard that splits any abelian 
subvariety up to isogeny) is given by the algebraic map to: 

p' 2 = 4p 3 - G 2 p - G 3 , G 2 = ^(g 3 2 + 9(7 3 2 ), G 3 = ^f-g 3 (3g 2 - g 3 2 ), 

P = -\ (4z 3 - 9g 2 z - 9g 3 ), p' = -- w (4z 2 - 3g 2 ). 

This map induces the following reduction of the holomorphic differential to an elliptic diffe- 
rential: 

dp = _^d_ 

p' w 

Notice then that in the equianharmonic case g 2 = the curve is singular, w 2 = — (4z 3 — 21g 3 )z 2 
and we have the following birational map to the elliptic curve: 

/2 , 3 ^ z , -2w „ 27 

p =4p -G 3 , P = -p P=^> Gs = -^ 3 

and differential dp/p' = zdz/2w. This motivated us to investigate the equianharmonic case 
further. 

What we believe to be new is the following observation. We found it surprising that there 
always be at least a second reduction in the case the elliptic curve is equianharmonic (the unique 
one that has an automorphism of order three), but see no theoretical reason to expect that. If 
this were to hold for all n, one could further ask whether these curves support KdV solutions 
elliptic both in the first and second time variables (a private communication by A. Veselov gives 
such an indication, to the best of our understanding). The issue of periodicity in the second KdV 
variable ('time') was put forth in |30| . and recently reprised in [10], but the double periodicity 
does not seem to have been addressed so far. We exhibit three covers in the g = 3 case (N.B. The 
third one is not, in general, over the equianharmonic curve); as well, we found two covers in 
the equianharmonic, n = 4, 5 cases, but we only provide a summary table which suggests the 
general pattern for two covers in the equianharmonic case. 

For n = 3, f s = z, fi = z 2 — Qzei + 45e? — I5g 2 , the curve is 



(c 



2 



(2376 z 3 g 3 - 36450 zg 2 g 3 + 504 g 2 z 4 - 91125 53 2 - 16 z 6 + 3375 g 2 3 - 4185 g 2 z 2 ). 



the first cover is given by the algebraic map 

(-84375c/ 3 ; + 2278125^1 + 3037505 3 ^ - 3375# 2 > 2 - 27000z 3 # 3 + 360z 4 g 2 + 16z e 



P 



36 z (-75g 2 + 4z 2 ) 2 
w (16z 6 - 1800,^2 - I6875gjz 2 + 421875s! + 54000z 3 # 3 - 11390625sf) 



108 z 2 (-75g 2 + 4z 2 Y 
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This map induces the following reduction of the holomorphic differential to an elliptic diffe- 
rential: 

dp of a 2 ir: \ dz 

— = -3(42 - 1552) — - 
d w 



In the equianharmonic case the above reduces to the following cover of the curve p 



,2 



4p 3 - g 3 - 



P 



1 (16z 6 - 27000z 3 £ 3 + 2278125c/!) 
~576 7 5 ' 

1 w(\Qz G + 54000z 3 # 3 - 11390625^; 
6912 ^ 



This map induces the following reduction of the holomorphic differential to an elliptic diffe- 
rential: 

dp = l2 z2 dz 



P' 



in 



Continuing in the equianharmonic case, the second cover is given by the following algebraic 
map to: 



P = 4p 3 - G 3 , G 3 



387420489 



P 



1 3645g| - 1080z 3 ff 3 + 16z 
6 z 



6 , _ 1 w(-729g 2 - 432z 3 g 3 + 16z e 

P 



32 



z< 



This map induces the following reduction of the holomorphic differential to an elliptic diffe- 
rential: 

dp = _ 10 ^£ 
p' w 

The third cover is given by the algebraic map to the (non-equianharmonic) curve 

p' 2 = 4 p 3 - G 2 p - G 3 = 4 (p - 198 53 ) (12717^ + 198p 53 + p 2 ), 

G 2 = 105948c/| , G 3 = 10071864^, p = 198^3 - 4z 3 , p' = 4wz. 

This map induces the following reduction of the holomorphic differential to an elliptic diffe- 
rential: 

cip zdz 
p' w 

For the curves: 



n 


fs 


Si 


1 


fs 


= 1 




fi 


= z - e. 


2 


fs 


= z 2 


- 3#2 


fi 


= z + 3ej 


3 


fs 


= Z 




fi 


= z 2 - 6zei + 45e 2 — 15g 2 


4 


fs 


= Z 3 


- 52g 2 z + 560^3 


fi 


= z 2 + Wza - 7g 2 - 35e 2 


5 


fs 


= z 2 


-27«? 2 


fi 


= z 3 - 15z 2 e 2 + (315e^ - 132g 2 )z + 675e| + 540^ 
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we summarize the equianharmonic covers in the table below (note: there is some small numerical 
discrepancy, e.g., for n = 2 our choice of normalization corresponds to the — 2w of this table, 
but we find it more useful to have the normalization below for recognizing a pattern). The pairs 
appearing in the second column are the orders of the commuting differential operators. 



n 


(",-) 


P 


P' 


dp 


G 2 


G 3 








P' 







2 (2,5) 9g 3 -Az 3 Awz 2 -3— 486#f 729^ 

dz 

3 (2,7) 198c/ 3 - Az 3 Awz -3z— 10594 5 f 10071864g§ 

dz 

4 (2,9) 1430^3 -Az 3 Aw ~ 3z2 — 6342300^ 3011499000#§ 

uo dz 

5 (2,11) 6435^3 -Az 3 4- -3z 3 — 28674000c/ 2 . 14791410000^ 



In factorized form the coefficients are 



(-,-) 


P 


P' 


dp 
P' 


G 2 


G 3 


(2,5) 


3 2 53- 


- Az 3 


Awz 2 


_ 3 dz 
w 


2 • 




3 6 5 3 3 






(2,7) 


2-3 2 


■ llg 3 - Az 3 


Awz 


dz 

— 6Z 

W 


2 2 


• 3 5 • 109p| 


2 3 • 3 6 • 


11 • 


157 5 3 


(2,9) 


2 • 5 • 


U-13g 3 -Az 3 


Aw 


-3z 2 ^ 

W 


2 2 


• 3 7 • 5 2 • 29g 2 3 


2 3 • 3 11 


•5 3 


• 17S 3 3 


(2,11) 


3 2 -5 


■ 11 • 13^ - Az 3 


A™ 

z 


, dz 
-3z 3 — 

W 


2 1 


■ 3 5 • 5 3 • 59g| 


2 4 • 3 6 • 


5 4 - 


2029^1 



The hyperelliptic curves are singular for g = 2, 5, 8; this is important because it makes it possible 
to explicitly calculate the motion of the poles of the KdV solutions (Calogero-Moser-Krichever 
system). When all the KdV or KP differentials (of increasing order of pole) are reducible, of 
which case we have examples in genus 3, both hyperelliptic as we have seen, and non (Section U]), 
the Calogero-Moser-Krichever system is periodic in all time directions; we believe this situation 
had not been previously detected. 



3.3 The theta function 

The next challenge we discuss is the calculation of the coefficients of the ODOs in terms of 
theta functions; while we do not present a result that could not have been derived by classical 
methods, we streamlined a two-step procedure: firstly, we use Martens' thorough calculation of 
the action of the symplectic group on a reducible period matrix [19] , to write theta as a sum, 
then we use transformation rules for the action under the symplectic group. The step which is 
not obvious however precedes all this and is the calculation of the explicit period matrix based 
on a reduced basis of holomorphic differentials. We give the result for one specific curve as an 
example. For the curve of genus 2: 

w 2 = (. 2 -e 1 2 )(. 2 -c 2 2 )(z 2 -ei), 

where we assume < £1 < £2 < £3, which is easily related to the form given by Jacobi for the 
case of reduction: 

w 2 = z(l — z)(z — a)(z — b)(z — ab), 
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our calculation shows that the associated theta function can be put in the form 




1 = I 63 \\ vu \ Ti ) °* (l vi - V2 ' \ T2 ) + \ 6i ll vi > \ Ti ) 6a \\ vi - W2 ' r 2 



where the Q{ are the standard g = 1 Jacobi theta functions, t\ is the r associated with the 
elliptic curve 

w \ = (x-e 1 ){x-e2){x-tD 

and T2 is the r associated with the elliptic curve 

w \ = x{x-i\){x-il){x-il). 

These two r's can be expressed explicitly in terms of elliptic integrals of the first kind K. 
We adopt the homology basis shown in Fig. 1. 




Figure 1. 



Ol 



We find that 

2^2,1 



2f2 2 ,i 



dx 

ui = I — , 



dz 

— = -w 2 , 
w 



zdz 



if 



20i, 2 

20 2 ,2 : 



-/ 


dz 

— = 1 


J 0.2 


w 


■J 


zdz 
= 


J a? 


w 


r dz 






= u' 2 , 


b 2 w 




-J 


zdz 


Jb 2 


w 



0, 



^ = 2f 

Ibl w 

where the uji are calculated from the the homology basis shown in Figs. 2 and 3, 

f& dx f$ dx . /*a dx 

1 - / , V 2 = , 1^2= ' 



U) 





Figure 2. 



Figure 3. 
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Proof 

Using the change of variable z 2 = 

dz _ r e2 2dz 
w ./_„„ w 



2Q 



1,1 



"i 



x, 2zdz = dx, we have 
dx 



(-:>, 



dx 



2dz 



a 2 



dz 
w 



f-i 



+ 



dx 



I y/x(x-&)(x-%)(x-%) 

dx 



7 

61 2dz 
w 



y/x(x-%){x-%)(x-€i) 



2n 



2,1 



"i 



zdz 
w 



~ e2 2zdz 

e 3 W 

dx 



2uj 2 , 

dx 



20 



2,2 



2zd,z 2zd.z 



"2 



o 



+ 



-e 2 )(*-e 2 2 x*-e 2 ) 

by symmetry under 



and 



2n' M 

20 12 

2n 2)1 

20 2 o = / 



dz 

i 10 
dz 
w 

zdz 
w 

zdz 
w 



(by symmetry), 



62 2dz 
w 



dx 



C 2 



2zdz 



w 



<-i 

62 2zdz 



•«' 



2 

2o»i (as for 0' 2 2 ), 



U2, 



The a; can be calculated by standard integrals 



i 



U>2 



I 

L0 2 





dx 




V(x~ 




CD 




dx 




V(x~ 


g)(s-£)(x- 
r/,r 


el) 


y/x(x 


-efx^-eix* 

dx 


-el) 




? v^-ei 2 )(x-ei)(x-ei) v(ei - e 2 )e 2 2 



(el 


-el)e 2 


(el 


-m 


(el 


-ef)el 


(el 


-e?)e 2 2 



so the period matrices we obtain can be written in the form 



n 



i 



-U) 2 ^2 



-ujL 
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so 



z I I 

U)\ U)y 

US\ 1 0J2 , 1 US 1 



2 us' 2 + 2 us' x 



2n 



1 1 



Following Martens, if we define 



m 



2 1 
10 



then 



US -1 



m 



2 [0 1] • [n n'] 



Now define 



J 



2 h 
-h 2 



T 









2 


1 


1 


-2 





2 





-1 





1 








1 






where T is symplectic, T ■ J ■ T f = J, and 



m • T 



10 
1-20 



is in (Martens) standard form, and 



m ■ J ■ m 







shows that the Hopf number is 2. 

The transformed matrix f is calculated by first calculating the 2x4 matrix 



[A B] = [I 2 rj-iJ-T)- 1 



then 



f = A~ l B 



~T~1 

2 

1 
2 



1 



22 + r 2 



where n = us[/usi. 

For this matrix, we use Martens' transformation formula, generalized to non-zero characte- 
ristics: 

1 

Vi Tl 







"0 0" 




r 


1 

[-2 °J 







-—V] 

2 



Vi - 2V2 

2(2 + r 2 ) 







"0" 




1- 


1 






.2- 



+ 9 



1 



Vl Tl 

2 ' 2 



8 



2 ' 2 



f 1 — 2f 2 2 
2 + r 2 '~2 + r 2 



e 



Vi — 2v2 2 
2 + t2 '~2 + t 2 



1 1 



1 



1 



1 



7^3 n T l #3 n^l ~ V 2 , ~T 2 + -#4 TT^li 77 T l ^4 7T W 1 ~ V 2, 77 r 2 



1 1 



1 



1 



in Jacobi notation, concluding the proof. 
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4 Genus three 

In [H], we calculated, using Halphen's ansatz as described in Section [T] (cf. [35] as well), the 
Halphen spectral curve, trigonal of genus 3, for the operator: L = d 3 — Ibpd — (15/2)p' = 0: 

, / 2 25 \ / 2 135 

We showed that this curve not only admits reduction, but also has Jacobian isogenous to the 
product of three (isomorphic) elliptic curves. Here we consider the covers in the slightly more 
general genus 3 curve w 3 = [z 2 — X 2 )(z 2 + A|), Ai,A2 £ M, to curves of genus 1 and 2 and 
derive the corresponding r matrices and reduction theory On occasion we indicate results for 
the special case \\ = (135/3) </ 3 , \\ = (25/A)g 3 : 

~ , / 2 25 \ / 2 135 
C 3 : w = U 2 + -^93) U 1 - — ff3 

This genus 3 curve covers the equianharmonic elliptic curve (p') 2 — 4p 3 — g 3 in three different 
ways and all entries to the period matrices are expressible in terms of the modulus of this curve. 
For the specific curve C 3 , we have three covers, all covering (p') 2 = 4p 3 — 33, as follows. 
The cover wi is given by 

_ w 2 (16z 2 + 8IOOC/3) j _ 2 z{l§z 4 - 19000z 2 c/3 - 759375#f ) 

P 



25 (4z 2 - I35g 3 ) 2 ' s 125 {\z 2 - 135# 3 ) 

with holomorphic differentials given by 

dz 3 dp 
w hp' 

The cover tx^ is given by 

\/5w , _ 4z 2 - 55g 3 

P "20^1' P " 80^g~ 3 

with holomorphic differentials given by 

zdz 3dp 

^ ~ iWy^p 7 ' 

The cover ^3 is given by 

1 w^/b{QAz Q - 80 g 3 z* - 53OO33V - 30375^1) 



P 



4800 (-g 3 ) 4 / s z 2 (4z 2 + 25g 3 ) 

-ipio 



5 1152000(4z 2 + 25g 3 )z 3 g 3 2 ' 

where 

pio = ( - 61509375 g 3 5 - 15187500z 2 5 3 4 - 700000z 4 5 3 3 - 240000z 6 5 3 2 
- 19200z 8 g 3 + 1024z 10 )\/l5 

with holomorphic differentials given by 

dz idpVl5 5 2 / 3 



w 



2 250p(-< ?3 ) 2 / 3 ' 
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The general curve C 3 is a covering of the equianharmonic elliptic curve given by the 
equation 

u 2 = + {X\ + X 2 2 ) 2 . 

The cover ir is given by 

= w, v = 2z 2 + X 2 , — X 2 , 

with holomorphic differentials given by 

2zdz d\i 
3w 2 v 

4.1 Reduction 

In this section we follow a similar approach to that of Matsumoto [20] in a genus 4 problem, 
who developed earlier results of Shiga [27] and Picard [22]. Write the equation of C in the form 

4 

w 3 = Y[(z - A fc ) = (z- Ai)0 + Ai)(z - iX 2 ){z + iX 2 ). 

k=l 

We fix the following lexicographical ordering of independent canonical holomorphic differentials 
of C3, du\ = dz/w, du2 = dz/w 2 , dus = zdz/w 2 , and will define the period matrix based on the 
branch cuts given in Figs. 4 and 5. 




Figure 4. 

Figure 5. 



4.2 Riemann period matrix 

We introduce the following vector notation: 

x = (x 1 ,X2,x 3 ) T = ( <j) dui,...,(b dui] . 
b = (61,62, b 3 ) T =1 f du 2 ,...,f du 2 ) , 



c = (c 1 ,c 2 ,c 3 ) T = I (h du 3 ,...,<b du 3 ) 

\J ai J C13 J 



T 
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Then the matrices of a and 6-periods read 

A = I i dui j = (x,b,c), (2) 

B=(/dtiij = ( P Hx,p 2 Hb,p 2 Hc), (3) 

where -ff = diag(l, 1,-1). The Riemann bilinear relation says 
x T ffb = x T Hc = 0. 

The Riemann period matrix r = AB~ l belongs to the Siegel upper half-space H 3 . 
By using the symmetries of the problem we can express all the x period integrals in terms of 
the two integrals (see Fig. 6) 

_ r Xi dz j _ r iM dz 

Jo w ' Jo w' 



Figure 6. Paths of / and J integrals on sheet 1. To get the corresponding integrals on the other sheets 
we multiply by p and p 2 respectively. The a and b contours can be constructed by pasting together 
various combinations of these integrals and their reflections in the real and imaginary axes. 



We have 



x x = (p 2 -!)(/- J), 



X2 



-Xl, 



X3 



-2J + p(J + 2p 2 I + p(J - I), 



where p = exp{2i-7r/3}. In the case when X 2 /X 2 = 5/27, the integrals simplify further, since 
/ = — J(l + 2p)/3. In this case (which we assume in all that follows) we have 



zi =2J(p + l), 
So we have 



[(ft') T tt T1 



X2 



-2 J 
-Pih 
■(P + I)ci 



-xi, x 3 = 2J(p-2). 

2J(3p + l) 



2J 

-piki 
-(p + l)c 2 



PlC3 



2Jpi -23 px 23{p - 2) 
bi b 2 63 

Cl C 2 C 3 
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where p\ = p + 1. We define 

r = ((n / ) T )" 1 n T - 

We find after much simplification, using the properties of p and the Riemann relations, that 



1 

79 



62p - 13 17 p + 13 -bp + 38 
17p + 13 62,0 - 13 5p-38 
-hp + 38 5/9 - 38 45/0 + 53 



This matrix [{Q') T £l T ] satisfies the reduction criteria as defined by Martens [19] . since if we 
define 



n = [2 J 2 Jp] , H = [1 0] , 



then 



i? • [(o') = n • m, 



where 



M 



-1111 
3 1 



We can transform M to standard form using the symplectic matrix 



S 



since 



M.S 



h- ' 


CO 


-1 







-6 


2 " 





-1 


-1 







1 


-1 





1 







-3 


-1 
















1 



















-1 





-1 





-2 







2 


3 


-1 




" -1 














" 







1 





-5 









Following Martens, the corresponding transformed r is 



T-S-J 



T 



1 1 

5 + 5^ 
1 
5 





1 

5 

1 1 

2 + 5 P 



2 P 



2 P+ 2\ 



Expanding a theta function defined with this f, again following Martens, we will get a sum of 
5 products of g = 1 theta functions with g = 2 theta functions. The g = 1 theta functions will 
have tau value r = (1/5) (1 + p) and the g = 2 theta functions will have a tau of 



"5 0" 




1 





1 1 

+ 

1 



2 + 5 P 



P 



1 

2 P 



2 P+ 2. 



5 
1 



25 

5p + y 



2 P 



5 
2 

3 1 

2 P+ 2. 
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Again following Martens, we can reduce this 2x2 matrix to standard form to get eventually 

- 11 1 1 - 

20 + 20 P ~4 

1 3 1* 
L "4 2 + l p \ 

Expanding each of the transformed genus 2 theta functions with this theta will give a product 
of genus 1 theta functions with r = (1/20) (11 + p) and genus 1 theta functions with r = 
16(3/2 + (l/4) p). So finally we have 5-4 = 20 terms, each containing a product of three genus 1 
theta functions (with fractional characteristics). 

5 Conclusions 

We contributed to the theory of spectral curves of ODOs with elliptic coefficients routine algo- 
rithms to calculate: 

• The algebraic equation of the curve (always, in principle); 

• The period matrix (only if the periods can be chosen suitably and there is an explicit 
solution to the action equations); 

• A reduction method for the theta function. 

What remains to be calculated (Challenge IV) is the dependence of the coefficients on 
the time parameters. This is more difficult because it involves expanding an entire basis of 
differentials of the first kind. In [9], we were able to find the time dependence by implementing 
Jacobi inversion, thanks to the Hamiltonian-system theory which describes the evolution of the 
poles of the KP solution [TT] . 
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